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Abstract. Molecules containing one or a few hydrogen atoms and a 
heavier atom (hydrides) have been predicted to trace FUV radiation. 
In some chemical models, FUV emission by the central object or pro- 
tostar of a star forming region greatly enhances some of the hydride 
abundances. Two massive regions, W3 IRS5 and AFGL 2591, have 
been observed in hydride lines by HIFI onboard the Herschel Space 
Observatory. We use published results as well as new observations of 
CH"'" towards W3 IRS5. Molecular column densities are derived from 
ground state absorption lines, radiative transfer modeling or rotational 
diagrams. Models assuming no internal FUV are compared with two- 
dimensional models including FUV irradiation of outflow walls. We 
confirm that the effect of FUV is clearly noticeable and greatly im- 
proves the fit. The most sensitive molecules to FUV irradiation are 
CH"*" and OH"*" , enhanced in abundance by many orders of magnitude. 
Modeling in addition also full line radiative transfer, Bruderer et al 
(2010b) achieve good agreement of a two-dimensional FUV model with 
observations of CH+ in AFGL 2591. It is concluded that CH"'" and 
OH^ are good FUV tracers in star-forming regions. 

1 Introduction 

, The formation process of high-mass stars is still disputed. In the scenario of 
continuous accretion, the young stellar object (YSO) starts to burn deuterium 
before reaching the main sequence. When this burning reaches the outer shells. 
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the YSO inflates and the surface temperature is reduced (Behrend & Maeder 2001). 
The same happens if the accretion rate exceeds 10~^Mq yr~^ and YSO radii may 
exceed 100 Rq (Hosokawa, Yorke, & Omukai 2010). A continuously accreting 
YSO moves in the Hcrtzsprung-Russcl diagram above and parallel to the main 
sequence. As the temperature reduction depends on the accretion rate, it depends 
on initial circumstances when the surface exceeds 10^ K, the temperature at which 
a significant fraction of the YSO luminosity is emitted in ionizing UV radiation. 
Alternatively, such radiation can be produced by accretion or jet driven shocks. 
X-ray emission of the YSO is limited to a small fraction of the Herschel beam for 
the distant massive objects studied here and is overwhelmed by far-UV (FUV, 6 - 
13.6 eV) if present (Bruderer et al. 2010a). 

We present results of the subprogram 'Radiation Diagnostics' of the Herschel 
Key Program 'Water in Star-forming regions with Herschel (WISH, van Dishoeck 
et al 2011). The goal of the subprogram is to explore the possibilities of identifying 
FUV and X-ray emission through chemistry in deeply embedded objects, where UV 
and X-rays cannot be observed directly dnc to a high attenuating cohimn density. 
Key questions are to quantify such high-energy radiation and their influence on 
star- and planet-forming regions. 

Here wc test molecular UV tracers in massive YSOs known to have Ultra 
Compact H II regions indicating that their UV emission has begun. We focus on 
hydrides, such as OH, CH, NH, SH, H2O and their ions, 0H+, CH+, NH+, SH+, 
H20^; and H3O"*". Many hydrides have a high activation energy in their formation 
reactions. However, high-energy photons such as FUV or X-rays heat and ionize 
the molecular gas. The hydride abundances are greatly enhanced by the high 
temperatures (Sternberg & Dalgarno 1995). Ionized hydrides are chemically active 
and can drive substantial chemical evolution. If their chemistry and excitation is 
understood better, they may become valuable tracers of warm and ionized gas and 
of the embedded phase in star and planet formation. 
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Table 1. Volume averaged (8000 AU diameter) fractional abundances relative to hydro- 
gen (Htot) derived from chemical model calculations. The no-FUV model is spherically 
symmetric, includes cosmic rays, but no internal FUV irradiation. The FUV model is 
two-dimensional and assumes protostcUar FUV irradiating the envelope including the 
walls of the outflows (from Bruderer et al. 2010a). 
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2 Model 

Molecular abundances were calculated in a time-dependent chemical network ir- 
radiated by FUV and X-ray radiation (Hollenbach & Tielens 1999; Stauber et al. 
2007; Bruderer et al. 2009). Here we use the template density models of AFGL 
2591 (van der Tak et al. 1999) and calculate the dust temperature from irradia- 
tion. We assume a chemical age of 4 x 10^ yr and a YSO FUV luminosity of 4x 10^^ 
erg corresponding to a bolometric luminosity of 2x10^* Lq and to an effective 
temperature of 3 x lO'' K. Abundances are determined in time and radial position. 
Finally the abundance is averaged over a volume weighted by density (Tab. [1]). 

Two methods to test theory with observations are used here. In a first and 
simple approach, the molecular abundances are estimated from observations and 
compared with models. In a second approach, the line emission of the molecules 
at the derived abundance are modeled by radiation transfer, convolved with the 
Herschel beam, and directly compared with observed line strengths. 

Two models are used: In the 'no-FUV model, only cosmic rays ionize the 
infalling envelope without additional high-energy irradiation by the YSO. It is a 
one-dimensional model with radial dependence on temperature and density, where 
the temperature model of Doty et al. (1999) is used. It assumes a cavity of 200 
AU radius in the interior. The 'FUV model includes FUV emission by the central 
object(s) and assumes irradiation of the walls carved out by the outflows. The 
modeled outflow has negligible column depth and the shape of the cavity allows 
direct irradiation of the walls. This irradiation both ionizes and heats a thin layer 
out to a distance beyond the Herschel beam on scales of 10000 AU. 
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Fig. 1. Line transitions of CH^ observed with Herschel/HIFI towards W3 IRS5. Zero 
LSR velocity corresponds to the LSR frequency of 835.1375 GHz (Miiller 2010). The 
systemic velocity of W3 IRS5, -38.4 km s~^, is shown by a vertical dashed line. The 
blue dashed curve indicates the contribution of the emission component. The red curve 
combines the contribution of the emission and absorption components fitted by the slab 
model. 
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3 Observations 

Hydride observations of two massive YSOs, W3 IRS5 and AFGL 2591, were re- 
ported by Benz et al. (2010) and Bruderer et al. (2010b), respectively. Figure 
[T] presents an additional observation for the first time. In both sources, the lines 
of CH+ have a P Cyg profile in the J = 1 — transition and an emission peak 
in the J = 2 — 1 transition at the systemic velocity. Combining the two lines, 
the abundance and temperature of the emission component, and the abundance 
of the ground state level of the absorbing component can be estimated using a 
slab model (Bruderer et al. 2010b). It is assumed that the absorption region is in 
front of the emission region, possibly associated with an outflow of the YSO. The 
best fitting model, indicated in Fig. [TJ yields a column density of 9.6 x 10^^ cm^^ 
at an excitation temperature of 43 K for the emission component, and 4.4 x 10^^ 
cm~^ for the absorption component referring only to the ground state. Note that 
in both W3 IRS5 and AFGL 2591 the CH+ J = 2 - 1 hue is relatively narrow 
(AV « 7km s^^) and symmetric unlike emission expected from shocked regions. 
Emissions and absorptions at Vlsr > —25 km s~^ may be caused by foreground 
clouds and are not modeled. 

4 Comparison of Observations and Theory 

4.1 Comparison of abundances 

Column densities N{x) are extracted by integrating line fluxes of a species x, 
neglecting re-emission or re-absorption of the final state. For H3O+ the total 
column density summed over all levels using the rotational diagram is derived 
(Benz et al. 2010). Only the line component corresponding to the systemic motion 
of the YSO and the emission component of CH+ are used in the comparison. 

In Fig. [2] observations are compared to model predictions of volume averaged 
abundances xi^) (Bruderer et al. 2010a). A dashed line N{x) — a xi^) (where 
a is considered a free parameter) is drawn such that the number of points above 
and below the line are equal (not counting upper limits above). In Fig. [5] (left), 
the column densities derived from observations are compared with the no-FUV 
model. If FUV is included in a spherical model, some molecules, such as CII+ and 
0H+ are enhanced an order of magnitude in abundance (Stauber et al. 2005). 
However, the average abundances of these molecules increase another five orders 
of magnitude or more if the FUV irradiates additionally the walls of the outflow 
cavity (Tab. [T]). The 2D model including FUV irradiation of the outflow walls is 
shown in Fig. [2] (right). It fits significantly better. 

4.2 Comparison of line intensities 

Bruderer et al. (2010b) have modeled the CH+ line emission of AFGL 2521 in 
the respective Herschelheam at 835 GHz and 1669 GHz (26.3" and 12.7"). FUV 
irradiation in a two-dimensional geometry enhances the strengths of the CH+ 
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Fig. 2. Comparison of column density derived from observations with model abundances 
calculated by Bruderer et al. (2010a) from a chemical model. A linear fit (see text) is 
indicated by the dashed line. Left: without internal FUV irradiation; right: internal 
FUV irradiation with a luminosity of 4x10^^ erg s~^ . 
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Table 2. Comparison of AFGL 2591 chemical and radiation transfer model calculations 
with Herschel/HIFI observations. The volume averaged model abundance refers to a 
radius of 20000 AU. The range of the FUV model indicates the range of free parameters 
(from Bruderer et al. 2010b) 

J = 1 — and J = 2 — 1 lines close to the observed values. The agreement is well 
within the range given by the uncertainties of the chemical reaction rates and of 
the model assumptions (Table [2]). 

5 Discussion and Conclusions 

Hydride observations of FUV emitting YSOs by HIFI onboard Herschel fit models 
including FUV irradiation considerably better than models without. Hence, FUV 
irradiation can be identified from molecular tracers. In star forming regions, CH"*" 
and 0H+ are the most FUV-sensitive species, and thus the best tracers. Therefore, 
the limitation to a single FUV tracer, such as CH+ or 0H+, may be sufficient. 
Here we have emphasized CH"*" which is more enhanced and easier to detect. 
CH+ emission probes the dense envelope; the order of magnitude fit of model 
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and observations strongly suggest that CH+ line emission is a method to trace 
enhanced FUV in star-forming regions. 

0H+ was observed mostly in absorption in W3 IRS5 and AFGL 2591, and 
only W3 IRS5 has a small P Cygni type emission. In Fig. [2] the 0H+ absorption 
component was used. A proper excitation analysis of that molecule would be 
needed for use as a reliable tracer. 0H+ refers to more diffuse and lower density 
gas. 

The other method compares abundances. Omitting effects of radiation trans- 
fer is clearly a coarse simplification. However, it includes contributions of various 
species that may originate from different locations and be better suited for quan- 
tifying the FUV flux in more complicated situations. 

We emphasize that the effect of FUV irradiation is strongly enhanced by a 
two-dimensional geometry increasing the irradiated surface area (Bruderer et al. 
2009). In these regions, the huge impact of FUV on the hydride chemistry in 
high-mass star forming regions is is not only caused by ionization, but also the 
result of gas heating. 
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